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Abstract 

We find that re-analyzing the LEP/SLC data with light superpartners and low 
a s (m|) ~ 0.112 yields a better fit to the data than the Standard Model, gives a sat- 
isfactory description of the Rb measurement, and gives a better fit to Alr. A large 
body of low energy (q 2 <§C m 2 z ) data and analyses provide compelling evidence for 
a s (m|) ~ 0.112. Global fits to LEP/SLC data in the Standard Model, however, con- 
verge on a value of a 8 {m 2 z ) ~ 0.126. Recently it has become increasingly clear that 
these should be viewed as incompatible rather than values that can be averaged. We 
investigate the possibility that new physics is causing the LEP high value. To this end 
we have conducted a global analysis of LEP/SLC data in the Standard Model and also 
in the Minimal Supersymmetric Standard Model. Several predictions could confirm (or 
rule out) the results of this paper: light chargino and stop, top decays into stop and 
neutralino, large Rb, large Alr, and a higher Jl%. We briefly discuss the implications 
of low a s for more fundamental high-scale supersymmetric theories. 
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Introduction 



Recently it has become increasingly likely that there exists a genuine and tantalizing discrep- 
ancy between low energy (q 2 -C M§) determinations of a s and the value of a s extracted 
from LEP/SLC data at the Z-peak. Shifman ||] has argued persuasively that the high 
value of a s (m|) ~ 0.126 obtained by fits to q 2 = m 2 z data is incompatible with the val- 
ues of a s (m 2 y) ~ 0.112 extracted from low energy observables and run up to the Z scale. 
Indeed, graphical demonstrations |2| of all the various determinations of a s clearly show 
an apparent systematic separation of a s (m z ) between the low energy data and the Z-peak 
data. 

In this letter we will assume as correct the plethora of extremely precise || low energy 
determinations of a s (m 2 z ) ~ 0.112. Then the extracted a s (m z ) from LEP/SLC must either 
settle to a lower central value with more statistics, or there is a systematic effect which causes 
LEP/SLC to fit to an inaccurately high value of a s (m z ). Our primary goal in this letter is 
to investigate whether a s {m 2 z ) extraction in a supersymmetric model can be substantially 
lower than the value of ct s (m 2 z ) determined from Standard Model fitting procedures, thus 
reconciling low energy and Z-peak determinations of a s (m z ). 

One way to think of this is as follows. The LEP/SLC data has been analyzed assuming 
the Standard Model is correct. If instead light superpartners exist, then a new analysis of 
the data is required. All output quantities will change. In particular, we find that a s (m z ) 
is allowed to decrease by about 0.01; Rb is now more consistent with the experimental data; 
agreement with A^r is better; and in general the global fit to the data is good. A number of 
other authors have also noted that if Rb is explained by new physics, then a s will decrease 
(See for example Refs. ||], ||). Before such an argument can be taken seriously, it is 
necessary to show that it is quantitatively large enough and also that it does not contradict 
other observables such as left-right asymmetries, forward-backward asymmetries and Mw- 
We have explicitly demonstrated these features. 

Gauge coupling unification and low a s 

Before continuing further, we should digress on a related question: Is a s (m z ) ~ 0.112 com- 
patible with simple grand unified theories? One of the early successes of supersymmetric 
grand unified theories was their ability to unify the gauge couplings (e.g., in SU(5)) and 
predict values of sin 2 9w an d a s{ m2 z) which were in accord with experiment. As the data 
and analyses got better, and the errors several times smaller, most upper limits on mea- 
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sured a s (m z ) started to drop. Simultaneously, super symmetry model builders refined their 
calculations and the theoretical lower limits on the predicted a s (m z ) rose. As it stands 
today, the lower limit on a s (m z ) is 0.126 in a simple SUSY GUT theory || (no GUT scale 
threshold effects, intermediate scales, or non-renormalizable operator effects) with common 
scalar and gaugino masses, and squarks bounded below 1 TeV. While this lower limit is 
compatible with the quoted (jj a s (m z ) from LEP/SLC data, it is not compatible with 
a s (m|) ~ 0.112. 

An a s (m z ) crisis is actually welcome because it demonstrates that we can learn about 
high-scale physics from weak-scale data. It leads us away from minimal models such as 
the CMSSM [|| which assume common scalar masses, common gaugino masses, and precise 
gauge coupling unification with a desert between the weak scale and the unification scale. 
This minimal constrained supersymmetric model cannot produce a s (m z ) below 0.126 or 
Rb above about 0.2168; it is a very predictive model. GUT scale threshold effects and non- 
renormalizable operators both modify || |k| simple notions of gauge coupling unification 
based on a continuous running of beta-functions from the low scale to the high scale, as 
do effects at intermediate scales that do not affect the perturbative unification OiJ. As 
low energy data gets better it starts to resolve gauge coupling palpitations near the uni- 



fication scale. Several authors || [TT| have used the lower a s (m z ) values to get insight 



into the form of possible supersymmetric GUT theories. This is in stark contrast to non- 
supersymmetric GUTs which have extreme difficulty rectifying the very large first-order 
problems of gauge coupling unification and proton decay constraints with second-order 
threshold corrections as well as keeping the weak scale and unification scale naturally 
separate. 

It has been suggested (l3|] that if one simply abandons the common gaugino mass as- 
sumption then low values of a s {m z ) can be obtained. While we fully agree with Ref. |l3|] 
on the importance of resolving this a s "crisis" , this is a dramatic approach, and a testable 
one. It is disquieting because in a simple GUT theory the gauginos must unify in a single 
adjoint representation of the GUT gauge group to preserve the gauge symmetry. If com- 
mon gaugino masses are discarded then gauge coupling unification also seems to be gone. 
In string theory, however, it is possible to have gauge coupling unification without having 
a grand unified group in four dimensions fll4j| . Usually it is assumed that the gauginos will 
unify as well but this is not necessarily required. What is required is the raising of the uni- 
fication scale from the typical scale of 10 16 GeV where simple SUSY theories want to unify, 



up to the string scale ~ 10 GeV. This is a non-trivial task 15], requiring the introduction 
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of additional states which affect the running of the gauge couplings. For these reasons, 
results based on simple GUT gauge coupling unification without gaugino mass unification 
are difficult to obtain in a theory. 

In this letter it is not our purpose to promote any specific notions of the GUT scale 
theory, and we do not attempt to provide any additional insight into how a more funda- 
mental high-scale SUSY theory could predict a low a s {m? z ). We shall focus instead on 
the low energy data, and demonstrate how fits to LEP/SLC Z-peak observables with light 
superpartners could give lower a s (m z ) than fits without superpartners. We know that by 
combining intermediate scales fllql , which do not hurt perturbative unification, with high 
scale threshold effects |], [ll|] we can construct a theory with the couplings and spectrum 
that we find in this work. 



Extracting a s in the Standard Model 

The values of a s (m z ) at the Z-peak are extracted, mainly, from two classes of observables: 
Thad and jet event shapes. The most important observables in the Thad 

class are Tz, 

-Riept = Thad/riept, and 0"h a d- The fits for a s (m z ) in the two approaches yield |7|, [If], 17 1, 

a s (m z ) = 0.126 ±0.005 from r\ a d observables, and 
a s (m z ) = 0.119 ± 0.006 from jet event shapes. 

The error in the a s (m z ) determination from T^ad observables is statistics limited. The 
error associated with the jet event shape measurements is mostly theoretical, since the 
non-perturbative effects of hadronization must be folded into the perturbative parton level 
jet correlations. Furthermore, the perturbative QCD calculations for the event shape mea- 
surements [18, 19] are not universally agreed upon, which compounds the uncertainty. We 



therefore cautiously ignore the jet event shape determination, which are in any case only 
la from the low values, and concentrate on the observables. 

In an effort to analyze all observables at LEP simultaneously in the Standard Model 
and in the minimal supersymmetric model we have implemented super symmetric loop cor- 
rections in Z0P0LE p3] and interfaced it with the CERN library minimizer MINUIT pi for 



a complete x 2 fitter. The observables that we use in our x 2 fit are Oj = T^, <7h a d> Rbi Rc, 
Alr, A b FB , Ap B , R\ ep t = Thad/Tiept, and A 1 ^. Next we fix the Higgs mass to a low value 
consistent with supersymmetry (rrih = 100 GeV), and let MINUIT find the minimum \ 2 f or 
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M t and a s {m z ). The x is defined as 



^theory _ ^expt^2 

X ^ y 



.2 



(AO 



cxptx 2 



All the values of 0* heory are calculated within a specific model and the better the match 
between theory and data the lower the x 2 • Using the Standard Model we find 

Aft = 167 ± 15 GeV 
a s {m\) = 0.123 ±0.005 

as the results of our x 2 fit to the observables. These results are consistent with the fits 
obtained by the LEP Electroweak Working Group @ corrected for a light Higgs. 

Extracting a lower a s in supersymmetry 

Now we set a s (m z ) to a smaller value (we choose 0.112) consistent with the numerous low 
energy observables, and map out the supersymmetric parameter space which yields a better 
X 2 with superpartners in loops and a s (m z ) = 0.112 fixed than does the Standard Model, 
whose x 2 minimum is at Mt = 167 GeV and a s (m z ) = 0.123. 

The idea that light superpartners might resolve the a s (m z ) discrepancy between high 
scale and low scale data is hinted at by the large measured value of Rb = T(Z —* bb)/T(Z — > 
had) which is approximately 2.3a from the Standard Model prediction. It was found in 
Ref. |^2[ that if and rn x + were both less than about 110 GeV then the discrepancy 
between theory and data for this one observable could go away. Since Rb had the high- 
est "pull" on the Standard Model x 2 f° r LEP data, resolving this 2. 3d deviation could 
substantially improve the global fit. 

If the theoretical prediction for Rb is raised by increasing the partial width, then for 
a fixed a s the total hadronic decay width is also increased. To a good approximation the 
hadronic width of the Z is separable into an electroweak piece and a QCD correction: 

■p theory _ -ptheory I , a s( m z) , | -pexpt 

L had — L EW,had I ~ l ~ n ^ J < ^ = ^ had 

Although Rb is rather insensitive to the QCD corrections, the partial widths Ti b and Th a d 
are quite sensitive. It is clear from the above equation that if we obtain a higher ^Ew^had 
in supersymmetry than was found in the Standard Model then the QCD corrections must 
be smaller in the supersymmetric theory to match the experimental determination of r^ad*! 
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that is, a s (m z ) must be lowered to best fit the data. Therefore, it qualitatively appears 
that we can simultaneously increase R^ and lower o. s , while at the same time keeping r^ a ^ ^ 
fixed. 

Our next step then is to hone in on the region of supersymmetric parameter space 



which will substantially increase Rb [22] and check to see that the x 2 fit to LEP/SLC data 
is consistent with low a s (m z ) and all other observables such as Apg, T z , R\epti etc. With 
light superpartners having a large effect on observables such as Rb, one would expect a 
priori that these same superpartners will affect other observables at LEP and potentially 
could yield a worse x 2 tit to the data than the Standard Model. It is imperative that 
all observables be analyzed simultaneously to confidently state that a lower a s extraction 
at LEP is possible in super symmetry. To be precise about our procedure, we have fixed 
a s (m z ) = 0.112 and searched through the MSSM parameter space for solutions which yield 
better x 2 , at fixed a s (m z ), than the lowest x 2 fit in the Standard Model where a s (m z ) was 
allowed to vary to its best-fit minimum value of 0.123. 

We have fixed a s (m 2 z ) = 0.112 for two reasons. One, we want to see if x'susy & t a 
low value of a s (m z ) ~ 0.11 can give a better x 2 than the Standard Model. And, we have 
determined that a s (m z ) = 0.112 is near the best minimum Xsusy ™ this analysis (with 
heavy first and second generation squarks and sleptons). Due to the extremely complicated 
minimization procedure with all the free MSSM parameters we do not yet claim with cer- 
tainty that the global minimum of the Xsusy fit is at a s (m z ) = 0.112, but only that there 
are at least local minima with a s (m z ) = 0.112 ± 0.004 and X'susy < Xsm- Furthermore, 
we have fixed tan /3 at its lowest possible value, which is determined by the top Yukawa 
remaining perturbative below the GUT scale, since this value gives the best Xsusy m the 
region of tan/3 < 30. For tan/3 > 30 the light pseudo-scalar Higgs can become important 
and we have not yet incorporated it into Z0P0LE. 

We have included into Z0P0LE all vector boson self-energy diagrams and vertex correc- 
tions which involve the charginos, neutralinos, stops and sbottoms. The only light squark 
or slepton expected in the spectrum which will affect our analysis is the tn, which becomes 
light through mixing in the stop mass matrix. Since the sbottoms are isospin partners to 
the stops they must be explicitly included in the calculation. We expect and assume that 
all other sparticles have masses too large to have a significant impact on the final answer. 
Although we work basically in a minimal supersymmetric theory, our results are largely 
independent of the gluino mass, and of first and second family squark masses if they are at 
all heavy. Results do assume Mi = M2 (bino and wino masses) at the GUT scale. Other 
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Figure 1: Region of supersymmetric parameter space with a better x 2 fit with a s {m 2 z ) = 0.112 than 
the best standard model \ 2 fit which was at a s (m%) = 0.123. 



Figure 2: Four observables versus the lightest chargino mass. The dotted line is the measured central 
value of the observable, and the dashed lines are the la limits. The solid straight line is the Standard 
Model best fit value obtained from Z0P0LE with rrih — 100 GeV, and the shaded region that which 
yields Xsusy < Xsm as other parameters are varied. As expected in a better \ 2 fit, the Rb and A^r 
predictions fit the experimental values as measured by LEP/SLC better than the SM does. Note 
also that the W mass prediction in supersymmetry is higher than the Standard Model prediction. 
And, the top is expected to decay into the lightest stop and light neutralinos with branching fraction 
as high as 60%. 

parameters are varied over allowed values (rather than guessed), to give the regions in the 
figures. 

Our calculations of the one-loop diagrams were checked in ZOPOLE by exact numerical 
cancellations of the log(// 2 ) which accompany all divergences in counter terms of the on-shell 
renormalization scheme. These exact cancellations of the log(fjp) in all observables and Ar 
are crucial requirements for a trustworthy calculation. 

Figure |l| is a summary of the main result in this letter. The enclosed area in the m x +~ 
mi plane is the region of parameter space which yields a better x 2 fit to LEP/SLC data 
using supersymmetry and a s (m z ) = 0.112 than the absolute lowest \ 2 obtained in the 
Standard Model (with a s (m 2 z ) = 0.123). The SUSY x 2 /d.o.f are as much as 1/3 better 
than the Standard Model best fit, and this minimum occurs when the chargino is near 
80 GeV and the stop is near 60 GeV. Interestingly, the lower bound on the lightest chargino 
is about 58 GeV although high i?^ values were obtained for m + < 58 GeV. The reasons 
for this are clear. The lightest neutralino in this region of parameter space is too light, 
and the Z decay width becomes too large. The truncated section in the lower right corner 
has a straightforward explanation as well. Here the stop is always lighter than the lightest 
neutralino and therefore becomes the LSP, which we exclude. 

It is very interesting to see the effect of supersymmetry on other observables. In Figure § 
we plot three observables, R b , M\y, and Alr versus the lightest chargino mass. The dotted 
line in each graph is the central measured value of each of these observables, and the dashed 
lines are the la errors associated with the measurements. The measured value for R b is 



taken from Ref. M, Myy from and Alr from [24|. The solid straight line is ZOPOLE's 
best fit Standard Model value with rrih fixed at 100 GeV (the Standard Model values would 
disagree more with experiment if rrih > 300 GeV). The shaded region is the range of values 



7 



obtained (versus lightest chargino mass), as other parameters vary, which yield a better x 2 
with light superpartners and a s {m 2 z ) = 0.112 than the best x 2 i n the Standard Model. 

Several aspects of Figure || are important. The Rb region is significantly higher than in 
the Standard Model. Myy is also higher. It is amusing that earlier values of M\y would have 
preferred the Standard Model to super symmetry, but the new value [^3| (80.33 ± 0.17 GeV) 
does not. The SUSY Alr value is closer to the SLC Alr measurement. These results 
translate to sin 2 0y/ = 0.2312 ± 0.0004. The values of M t that we found with x'susy < Xsm 
range between 162 GeV and 190 GeV. The upper limit on Mt comes about mostly from 
the inability to get low tan (3 and high Mt simultaneously, and still keep the top Yukawa 
perturbative at the high scale. With very light charginos we run the risk of having top 
decays into the lightest stop and light neutralinos be too numerous to be consistent with 



top quark production and decay data at Fermilab [25]. Figure || shows that the branching 
fraction of these supersymmetric top decays can be as high as 60%, and in general much 
of the parameter space has a significant top decay branching fraction into supersymmetric 
states which could be detected when many more top events are detected at a high luminosity 
collider. 

It should be re-emphasized that the most important phenomenological implication of 
lowering the extracted a s {m z ) is light superpartners. Most of the allowed parameter space 
in Figure 1 will be detectable at LEP II and an upgraded FNAL collider. With sufficient 
luminosity LEP II will be able to detect all charginos and stops with masses to within a few 
GeV of yfs/2. An upgraded Tevatron collider should be able to reach charginos and stops 



with considerably higher masses [26] than LEP. However, FNAL, and to a limited extent 
LEP, has some difficulty cleanly detecting a signal for Higgsino-like charginos. In the limit 
of pure Higgsino the LSP mass gets closer and closer to the lightest chargino mass. When 
the chargino decays into LSP plus leptons, the leptons may have too little energy to trigger 
on, so the signal is reduced. This region of chargino parameter space is largely the region 
we are in. 



Conclusion 

We have demonstrated that the extracted value of a s (m z ) from LEP/SLC data can be 
lowered to agree with other a s (m z ) determinations when superpartners are added to the 
fit. An essential aspect of this work is the inclusion of all relevant LEP/SLC data, so 
that the results are known to be consistent with all observables. We have found that light 
charginos and stops (with masses below ~ 100 GeV) are required if the total x'susy with 
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a s {m 2 z ) = 0.112 is better than the Xsm with a s( m2 z) a t its Standard Model best-fit value 
of 0.123. Our approach is largely independent of SUSY assumptions. 

The SUSY spectrum and couplings required to obtain our results cannot be obtained in 
a fully minimal supersymmetric model. They can be obtained by adding the effects of high 
scale thresholds, and/or Planck scale operators, and/or perturbatively valid intermediate 
scales. It is very encouraging that data at the electroweak scale seems to be telling us about 
physics near the Planck scale. 

The resultant supersymmetry parameter space has several important phenomenological 
implications: The W mass is higher than the expected Standard Model best fit. Rf, and 
Alji should also be larger than their Standard Model values. Light superpartners below 
about 100 GeV must exist. LEP II and FNAL will probably find these superpartners if they 
are this light; if they don't, very precise determinations of the W mass, Rb, or A^r could 
rule out or further support this exciting possibility. 
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